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The formation of stable dimers on top of terraces during epitaxial growth is investigated in detail. 
In this paper we focus on mean-field theory, the standard approach to study nucleation. Such 
. theory is shown to be unsuitable for the present problem, because it is equivalent to considering 

adatoms as independent diffusing particles. This leads to an overestimate of the correct nucleation 
rate by a factor A/", which has a direct physical meaning: in average, a visited lattice site is visited 
Af times by a diffusing adatom. The dependence of N on the size of the terrace and on the strength 
of step-edge barriers is derived from well known results for random walks. The spatial distribution 
of nucleation events is shown to be different from the mean-field prediction, for the same physical 
reason. In the following paper we develop an exact treatment of the problem. 
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I. INTRODUCTION 

A crystal can be produced artificially with different techniques: we can pull it from the melt, grow it from a solution 
or obtain it via deposition from a gas/vapor phase onto a suitable substrate. One of the key mechanisms of the growth 
process is the formation of supercritical nuclei, that is the nucleation - via diffusion and aggregation - of crystalline 
\ clusters whose growth rate exceeds the decay rate. 

In this paper we devote our attention to epitaxial growth by atomic or molecular beams particles travel 
ballistically towards the growing surface where they undergo a thermally activated diffusion process. The size i* of 
the critical nucleus is typically a few units and its actual value depends on several factors: the substrate and the type 
of adatoms deposited determine the activation barriers for the different atomistic processes, while the temperature 
7-H , and the flux determine what processes are really relevant on the time scale of the experiment. Here we consider 
the simplest case: nucleation is irreversible (i* = 1), i. e. once two adatoms meet they form a stable dimer. After 
j nucleation, the stable nucleus grows by capturing other adatoms. 

We start by explaining in general terms the role of nucleation in the different stages of epitaxial growth and by 
discussing how rate equations and mean- field theory deal with it. The focus of the rest of the paper will be on 
. nucleation on top of terraces bound by descending steps (often called second layer nucleation) . 
' In the submonolayer regime diffusion takes place on the substrate: adatoms are deposited randomly and they diffuse 
f^) , until they meet another wandering adatom or a growing cluster. 

Rate equations |2j are widely used to describe the processes of adatom capture and adatom detachment from clusters 
of j atoms (j— clusters, j > 2). If only adatoms are mobile, the spatial density pj of j— clusters varies in time according 
to the relation [|| d t pj — Uj—i — Uj where Uj is the net rate of the process (j— cluster)^ (j + 1-cluster). Uj is the sum 
of a "growth term" and a "decay term" : the growth term represents the aggregation of an atom into a j— cluster and 
it has the form ajDppj, where p is the adatom density, D the adatom diffusion constant and Oj is an adimensional 
capture number; the decay term represents the thermal detachment of an atom from a (j + 1)— cluster and it has the 

form — (^fj Pj+i °xp(— Ai?/T), where ao is the lattice constant and AE is the energy difference between a j— cluster 

(plus a free adatom) and a j + 1-clustcr. If i* = 1 such term is absent, because all j— clusters are stable and the 
nucleation rate is 07 RE = <J\Dp 2 (the subscript standing for "Rate Equations"). 

The capture factor Uj is defined through the flux $>j of adatoms attaching to the j— cluster: $j = Dajp, and it 
accounts for the different adatom densities surrounding islands of different size. In a Mean-Field (MF) approach such 
density is taken as a constant, a — 2d H and the nucleation rate per lattice site reads uj m¥ = 2dDp 2 . 

By monitoring the adatom density p and the total density of islands p TOT = X)j>2 Pj one realizes the existence of 
a regime characterized by p TOT almost constant in time, prior to the coalescence regime [|J. The average distance t-o 
between islands can therefore be defined: £^ = 1//0 TOT , d being the dimensionality of the substrate: usually d = 2, 
but in the following we will consider d = 1 as well. £ D is called the diffusion length and it also gives the typical linear 
distance travelled by an adatom before being incorporated in an island. 
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During the time 1/F, necessary for the deposition of a monolayer, there is on average a nucleati on e vent per area 
i v |6fl. Therefore, according to mean-field theory, lu mf ■ i*/F ~ 1, i. e. DflijF ~ 1 with (see App. [All) p ~ So 
we obtain the result £ D ~ [D/F) 1 , with 7 = ^ in d = 1 and 7 = i in d = 2. This relation is wrong in one dimension, 
but essentially correct in two dimensions In Sec. IV A it will become clear that Mean-Field Theory (MFT) gives 

the correct result for 7 in d — 2 because step-edge barriers at island edges play no role in the submonolayer regime. 

If adatoms freshly attached to a growing nucleus are immobile the resulting island is normally fractal, but if they 
can diffuse along the step-edge the island gets compact. Compact islands are therefore obtained at not too low 
temperatures (so that thermally activated edge diffusion actually occurs) whilst nucleation is irreversible at not too 
high temperatures (so that a dimer is thermally stable): in some experimental systems both conditions are fulfilled 
in a given temperature range: Pt(100]/Pt in the range 350-430 K g, Ag(100)/Ag in the interval 200-300 K and 
Fe(100)/Fe between 300 and 500 K fnj] are some examples. 

Once island coalescence has set off, most nucleation events take place on top terraces. The obvious reason is 
that even for layer-by-layer growth a new atomic layer starts to form before the completion of the previous one and 
consequently the growing surface is made up of a certain number of exposed layers. Terraces can be classified as 
top (T), vicinal (V) and bottom (B) terraces according to the type of steps surrounding them: in general only top 
terraces attain a size large enough to have a considerable probability of nucleating a new island. Since this probability 
grows abruptly from zero to one with increasing terrace size , it is possible to introduce a critical nucleation length 
Cc |J- In the presence of step-edge barriers |13| , hindering interlayer transport, such length differs for the three types 
of terraces: £^ uc goes to zero with increasing barriers while £^ uc and £® uc remain finite in such limit |T^,Q. Nucleations 
generally occur on top terraces because the size of vicinal and bottom terraces hardly reaches the nucleation length. 
However nucleation is a stochastic phenomenon and therefore dimer formation may also occur on vicinal terraces every 
now and then, while on bottom terraces it is an exceedingly rare event. These occasional "vicinal nucleations" are an 
important stabilizing mechanism of the surface, even if their relative weight decreases as barriers increase 

The mean-field approach to evaluate the nucleation rate for (reversible and irreversible) nucleation on top terraces 
has been worked out by Tersoff, van der Gon and Tromp in Ref. [jlTf and since then their results have been widely used 
to analyze experimental data and extract the values of step-edge barriers. The mean-field approximation can also be 
used in 'mesoscopic' models of growth fl7| for implementing the rule for the spatial distribution P(x) of nucleation 
events, via the relation P(x) cx p 2 (x). 

Recently, several authors |T^-|2b]| have started to review critically the MF approach. In Refs. []l8|,|o| authors are 
mainly interested in the validity of MFT for different values of the critical size i*: this check has been done through 
scaling analysis and kinetic Monte Carlo simulations. A deeper investigation of the irreversible case (i* = 1) is done 
in Ref. [ 1^] in the limit of strong step-edge barriers, through a quantitative approach based on the different time scales 
involved in the nucleation process. 

The inaccuracy of MF theory in dealing with the problem of second layer nucleation has therefore already been 
pointed out in the literature. In this work we analyze thoroughly the validity of MF theory for irreversible second 
layer nucleation, we rigorously prove its inadequacy both for the nucleation rate lo and the spatial distribution P(x) 
of nucleation sites and understand the physical origin of its failure. In the following paper pl| we compute exactly lo 
and P(x), providing the correct expressions that must be used in place of the MF approximations. 

We have organized this paper as follows. The three relevant time scales for the nucleation process and the associated 
regimes are introduced and discussed in the next Section. In Sec. [II we explain how we can get rid of the stochastic 
nature of the depositi on pro cess (III A) and we introd uce th e quantities of interest in the paper, the spatial distribution 
of nucleation events ( [II B ) and the nucleation rate ( III C ). In Sec. [V we show the equivalence of MF theor y with 
a model of noninteracting particles, both for the nucleation rate (IV A) and for the spatial distribution (IV B). The 
final Section contains a critical summary of the results. 

A short report of this work has been published in Ref. G2J . 



II. TIME SCALES 



We consider a top terrace of fixed linear size L, subject to a flux F of particles. Once on the terrace, each particle 
moves with a diffusion constant D until it meets another particle or it leaves the terrace. At step edges an additional 
energy barrier, usually called Ehrlich-Schwoebel (ES) barrier p3| , reduces the rate D' of interlayer transport and the 
ES length 

% i) «o (1) 
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measures the asymmetry between D and D'. In the following the lattice constant ao will be taken as unit length 
(a = l). 

In general, throughout the paper we will consider discrete space and time, i. e. particles moving on a lattice (a 
square lattice in d = 2), at fixed time steps. However we will sometimes use a continuum notation as well. The 
matching between discrete and continuum is straightforward: if t c is a continuous time and td its discrete counterpart 
(number of time steps), then t c = tj,- At, where At = 1/(2gLD) is the time an atom remains on each lattice site. 

Let us now discuss the three time scales involved in the problem Jl9[] . 

• The traversal time (r tr .) is the average time needed by a diffusing particle to reach the terrace boundaries. In the 
large L limit 

T tr = PL 2 /D , (2) 

where (3 is a numerical prefactor depending on the dimension d and the shape of the terrace. Its value is discussed 
after Eq. (@). 

• The residence time (r res ) is the average time a particle spends on the terrace. It is related to the average density p 
of adatom s vi a the relation [l9| ] p — Fr res , where the density p (and therefore its average value p) can be determined 



(see App. AT) by solving the stationary diffusion equation in the presence of a constant flux F: 

DV 2 p + F = 0. (3) 

Boundary conditions depend on the strength of the Ehrlich-Schwoebel (ES) barrier at descending steps. If the 
equilibrium adatom density - due to thermal detachment from steps - is neglected, the boundary condition for p is: 
d±p\ st op = j^rUcp, where d± is the local derivative in the direction perpendicular to the step (directed inward the 
terrace). In a discrete picture (for example in d = 1) if n = 1 is a lattice site at the edge of the terrace and n = is 
its fictitious neighbour outside the terrace, we have p(l) — p(0) = p(0)/£ ES , i. e. p(0) = 1 ^| s p(l) = ap(l). 

The solution of Eq. (||) (see App. Al) gives the following result, valid both in one and two dimensions: p = 



j) (PL + a£ ES )L, where a is another numerical factor depending on the dimension d and on the shape of the terrace. 
We can finally write: 



{PL + a£ ES )L/D . (4) 



In the absence of ES barriers (£ ES = 0) r res and r tr are equal. From Eqs. (|A7|,[A^, |Al8[|A5lD we infer that in d = 1 



P = 1/12 and a = 1/2. In d = 2, for a circular terrace of radius L, p = 1/8 and a = 1/2, while for a square terrace, 
P ~ 32/tt 6 and a = 1/4. 

In the discrete picture r res is clearly equal to the total number of sites (N all ) visited by an atom during its diffusional 
motion on the terrace. Since the adatom stays on a lattice site a time l/(2cLD) we have T res = N all / (2dD). Hence, the 
residence time r res , the average density p and the number N aU of all sites visited by the adatom carry the same piece 
of physical information, once F and D are set. The quantity A^ u should not be confused with the number of distinct 
sites {N dis ) visited by an adatom: a given lattice site, visited k times, contributes for 1 to 7V dis and for k to iV aU . 
• The deposition time (rdep) is the average time between a deposition event and the next one. For a terrace of area 
A = L d : 

11 

Tdep ~ FA ~ FIA ■ (5) 

Physically sensible values for F, L and D imply that r tr -C Td ep - This relation indeed is <C jqj, i- e. y ^> L d+2 : 
we can now recall the diffusion length introduced in Sec. |, £ D ~ {y) 1 an d measuring the 'maximal' size of a 
terrace in the absence of step-edge barriers. For irreversible nucleation the exponent 7 is equal to |23]| 2 (d+i) so that 

2(d+l) 

we obtain the condition £- D d+2 L, i. c. £^ 3 L in d — 1 and t^J" 2 L in d — 2. Smooth growth requires that 
£ D 3> 1; furthermore L is at most of order £ D if £ ES — 0, but for finite barriers it is (much) smaller. We conclude that 
the above conditions are fulfilled and that we can safely suppose that r tr <C Td ep - 

This inequality has a consequence of primary importance: processes involving more than two adatoms at a time 
can be fully neglected. Once two adatoms are simultaneously on the terrace, they meet - if they do - on the time 
scale of the traversal time r tr . This fact is intuitively clear and it is proven in Ref . ]2l| ] . The probability that a third 
atom lands in the meanwhile is T tr /Td e p, negligibly small. Irreversible nucleation is therefore the result of two-adatoms 
processes only. 

Depending on the relative size of r res with respect to the other two time scales, three different regimes may occur: 
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i) T tr — T res <C Tdep Zero or weak barriers (6) 
it) T tr "C T T es "C T<j e p Strong barriers (7) 
Hi) T tr "C T^ep "C T res Infinite barriers (8) 

The difference between the three regimes is easily understood. A nucleation may occur only if a new adatom is 
deposited before the previous one leaves the terrace. If T res <C Tdep (regimes i and it) this is a rare event. When it 
happens, the second atom finds the first one with a spatial distribution that differs in cases i and it (see App. |§|). If 
£ ES <C L (regime i) when the adatom reaches the edge of the terrace, it gets off. Steps act as absorbing boundaries 
and the adatom density vanishes there: p has a parabolic shape with a maximum in the middle of the terrace (see 
App. |A 1| ) . If £ ES L (regime it) the adatom is pushed back from the terrace edge several times before being able to 
descend. Steps act as (imperfect) reflecting walls and p is approximately uniform over the whole terrace. In regime Hi 
when a new adatom is deposited it always finds the previous one still on the terrace, they both have a flat distribution 
and they will certainly meet. 



III. GENERAL FORMALISM 



In the present and in the following paper we are going to use a discrete formulation for particle dynamics, both in 
space and in time. In this Section n indicates the whole set of d integer numbers specifying the position of a particle 
on the terrace. A nucleation event is assumed to occur when two adatoms are on the same lattice site, rather than 
on neighbouring lattice sites: this definition avoids useless mathematical complications, but retains all the physics of 
the nucleation process. 



A. Reduction to two particles deposited simultaneously 

It is clear that the problem of dimer formation on a terrace involves the study of the diffusion of two particles 
deposited at different times: the spatial and temporal distributions of landing events play therefore a prominent role. 
The incoming flux of particles is supposed to be spatially and temporally uniform [j24| : a particle arrives on each 
lattice site with uniform probability p^ = 1/L d and the interarrival time r between two deposition events decays 
exponentially p5[| : 

-Pdcp(T) = exp{-T / Tdep) / Tdep ■ (9) 

Let us now consider any "two-particles" quantity O, i. e. any quantity depending on the initial distributions of 
particle (1) and (2) and on their interarrival time r. 

Let p n {0) = Pn — ^/L d be the initial uniform distribution of an atom and p n (r) its dynamical evolution at time 
r (in the absence of other particles). If particle (1) is deposited at time zero and particle (2) a time r later, we call 
O{pJ l 1) (r),p<j 2) (0)} the resulting physical quantity. O might be, for example, the probability p nuc that a deposited 
particle nucleates a dimer before getting off the terrace [see Eq. (jl^)]. Once O{p^ 1) (r),p[ l 2) (0)} is known, one should 
evaluate its average over r: 

oo 

6 = Y J P^{r)0{p^(r)^ ) m ■ (10) 

T=0 

The crucial point is that if O is linear in the initial distributions p£' 2) of the two atoms (as all quantities discussed 
in the paper are), the previous equation can be rewritten as 

6 = o{pf,p^m (ii) 

in such a way that the average over r is now included in an effective initial distribution: 

oo 

iC = E p -pWM • (12) 

T=0 

We can make more explicit the physical content of the above reasoning, that is based on the linearity with respect 
to the initial distribution of the two particles. Atom (2) is deposited with probability P dop (T) a time r after atom (1), 
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which means that atom (2) has the probability P dep (r) to find atom (1) distributed according to P^'(t): on average - 
and it's now that linearity comes into play - atom (2) finds atom (1) with the effective distribution given in Eq. (|12j). 

In this way we have reduced the problem of evaluating O to the evaluation of O for two particles deposited 
simultaneously. Thus we can ignore the stochasticity of the deposition process and assume that atoms (1) and (2) 
land at the same time, but the actual initial distribution for atom (1) (the uniform distribution) is replaced by p^ ff . 

The next task is then the determination of p°". The function p n {r) (discussed in App. [A 2| ) is the distribution of 
the first adatom at time r, i. e., the solution of the diffusion equation for a single particle with initial condition 

p„(r = 0)=K = ^j • (13) 
The sum of p n (t) ove r all times r is the solution of the stationary diffusion equation ([?]) (of its discretized version, 



actually, see App. A 1), whose normalized form will be indicated with p s n (s standing for stationary). It has in general 
a parabolic form, and in particular in d = 1, p„ = [£ ES L + (L + l)n — n 2 ]/[£ ES L 2 + \L(L + 1)(L + 2)]. 
In Appendix [b| it is shown that in all dimensions we can write 

Pn = ~ ~~rZ Pn ■ (14) 
'dep t 'res 

The physical content of Eq. ( pT| ) is readily understood. For infinite barriers (regime Hi) p^ 1 — = 1/ ' L d : the first 
particle cannot escape from the terrace and its distribution is still uniform and normalized when the second one lands. 
For strong but finite barriers (regime it) p"" = ^p^-p^ = ^F^- : most of the particles that arrive on the terrace leave 
it before another particle lands, but the distribution of the first particle remains practically uniform because many 
attempts are needed to overcome the ES barrier. In the limit of zero or weak barriers (regime i) pff = ^-^pf, and pf 
vanishes on the edges, reflecting the presence of the absorbing boundaries. 



B. The spatial distribution of nucleation events 



In the previous subsection we have explained how to transform the original problem into the new problem of two 
atoms deposited at the same time, with normalized distributions p^ (the first) and p^ (the second). We can now 
define the probability R(n, t) that a nucleation event occurs on site n at time t and introduce the following quantities: 

P(n)=Y,R(n,t) , (15) 
t 

W = ^P{n) . (16) 

n 

P(n) is the spatial distribution of nucleation events and W is the probability that two atoms, both on the terrace 
at time zero, meet before leaving the terrace. It is useful to consider the normalized spatial distribution P (N) (n) = 
P(n)/W as well. 



C. The nucleation rate 



The nucleation rate lj is defined as the number of nucleation events per unit time on the whole terrace of size L, 
irrespectively of the spatial location of the meeting point. This quantity is of great importance because it is related 
to the probability of second layer nucleation. In a classical experiment p6]| a fraction of a monolayer is deposited 
on the substrate and the size of islands is made as uniform as possible through an annealing procedure. Starting 
from this template a second dose of atoms is deposited and nucleation on top of existing islands is monitored. oj(L) 
enters in the interpretation of this experiment because the probability P(t) that a nucleation event has occurred on 
a terrace by time t is V(t) = 1 — exp{— J Q drw[L(r)]}. The rate io(L) is defined and evaluated for a constant terrace 
size L: in the experiment instead, L grows in time and the time dependence of L is 'system-dependent'. Hence, the 
growth law L(t) of the terrace size must be supplied beyond to(L), and it depends on the specific morphology of the 
surface and the experimental setup. In other words, the nucleation rate - on the one hand - has a very general and 
basic meaning, but - on the other hand - it can hardly be measured directly. This means that, despite our results for 
to(L) are exact, the evaluation of L(t) introduces some approximations in the interpretation of experimental results, 
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whose accuracy depends on the detail of the system considered. In addition, some secondary effects, as steering and 
nonuniform barriers, may further complicate the problem. 

We now connect to to P(n) and W. In Sec. || we explained why only processes involving two adatoms are relevant 
for studying irreversible nucleation, because r tr <C Tdep- If we define the nucleation probability per atom, ]3 nuc , we 
can write the nucleation rate as the number of atoms landing on the terrace per unit time (FL d = l/r^p) times the 
nucleation probability per atom: 

uj = FL d -p mc . (17) 

The quantity p nuc is the probability that a deposited particle nucleates a dimer before getting off the terrace and it 
can be written as 

oo 

P„„ c = J2 P ^P™iPn^Pn(°)} > ( 18 ) 

T = 

where the dependence on the initial distributions of atoms (1) and (2) has been made explicit. We stress that the 
dependence on the initial distributions occurs via the full diffusion process. For example, for independent particles in 
one dimension, the explicit form of p nuc is 

L oo 

p nuc K l5 (T),K 2 '(0)} = £ ^ Pra , m (f)K'(r),^'(0)} (19) 

m=l t=0 

where p m ,m(t) is the solution of the diffusion equation in two dimensions with initial condition given by the product 
p(i)( r )p(2)(Q) ( S ee Ref. |lj] for more details). 
Because of the linearity of p nuc , we have 

p nuc = P „ u jK'V< 2, (0)} • (20) 

The nucleation probability per atom can be thought as the probability that atom (1) is still on the terrace when 
atom (2) is deposited, times the probability they meet before getting off the terrace. This is exactly what emerges 
from Eq. J20| ) once expression (|l4) for the effective distribution p°" is inserted: 

•p nuc K, P < 2 »(0)}. (21) 



Tdep H~ T r 



The normalization factor of p£f is the probability that atom (1) is still on the terrace when the next one shows up; 
for infinite ES barriers (r res 3> Tdep) such a probability is trivially 1, while for weak and strong barriers (r res <C Tdep) 
it is T res /Tde P - The remaining quantity on the right hand side (p m <>{Pn, PnK®)}) ^ s the probability that two atoms, 
both on the terrace at time zero [p® and p^ (0) are normalized] , meet before descending. Therefore it coincides with 
W and we finally obtain 

u = FL d ^ W . (22) 

Tdep T" T re s 



D. Noninteracting particles 

We are considering a system such that once adatoms come together an immobile dimer is formed irreversibly: 
adatoms stop diffusing and the dimer does not dissociate. It turns out to be of great help to consider also an artificial 
model, with adatoms treated as independently diffusing particles: even if they meet on the same lattice site they go 
on diffusing and therefore they can cross each other several times before leaving the terrace. We consider all these 
meetings as 'fictitious nucleations', and define also for noninteracting particles the quantities mentioned above: the 
nucleation rate w NI , the spatial distribution P m (n) and the total number W NI of nucleation events, the subscript NI 
standing for "Non Interacting" . 
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FIG. 1. (a) Log-log plot of W (empty circles), N diB /L 2 (full circles), Wm (empty squares), N aU /L 2 (full squares) versus £es 
for d = 2 and L = 20. (b) The same quantities plotted versus L for d = 2 and 1-es = 0. 



IV. EQUIVALENCE OF MEAN-FIELD THEORY AND NONINTERACTING PARTICLES MODEL 



A. The nucleation rate 



We have introduced the nucleation rate in Sec. Ill C and obtained Eq. (|22j). W is the nucleation probability between 
two atoms that are both on the terrace at time zero. For the model of noninteracting adatoms W should be replaced 
by Wnij the average number of meetings between the two independent particles. Of course Wm can be larger than 1. 

The simplest and less interesting case is the regime Hi of infinite barriers. In such a case W is trivially 1 and 
u> = FL d = 1/rdep, i- e. any particle deposited on a terrace does form a dimer. In a sense, this limit is unphysical for 
mean- field theory because p and u> MF — 2dDL d p 2 diverge when £ ES — * oo. 

In the other two regimes (weak and strong barriers), Eq. (22) becomes 



u = FL d '^L W = FL 2d pW , 

7~dep 



(23) 



where we have used the relations Tdep 
particles and obtain 



(FL d ) 1 and Fr res = p. We can repeat the same procedure for noninteracting 



w m = FL d — W m = FL 2d pW N 

Tdep 



(24) 



It is possible to relate W and Wm to single-particle quantities, N diB and iV all (see Sec. ||). They are the number of 
distinct (iV dis ) and all (iV all ) sites visited by a single walker diffusing on the terrace p7| . Let us assume one of the two 
adatoms fixed on site s. w(s) is the probability that the diffusing adatom visits site s before getting off the terrace. 
W is then the average value of w(s), W ~ J2 s w ( s )/L d - The quantity J2 S W ( S ) ^ s nothing but the total number 
of distinct sites N die visited by the diffusing adatom, so that W ~ N dis /L d . The same argument for noninteracting 
particles gives W N1 ~ N all /L d . 

The relations: 



W 



N dls 
L d 



and W N 



L d 



(25) 



have been derived under the assumption that one atom is immobile. In Fig. [I] we compare numerically the values 
of single-particles quantities (^S : ^ 1 ) with two-particles quantities (W, W N1 ). The former have been calculated via 
Monte Carlo simulations and the latter through the numerical solution of the discrete diffusion equation for two 
atoms on a terrace (discussed in detail in the next article): it comes out that relations (^5j) are well satisfied, so that 
assuming one atom as immobile is perfectl y r easonable for the evaluation of W and W m . 
If we insert the relations (p5|) into Eqs. (|23|,p4]) we obtain: 
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FIG. 2. Plot of the correction factor ui MF /u> and of Af = A/ln/A/dis versus £ E s/L, for a square terrace of size L = 20. 



FL pN d 



w NI ~ FL pNz.ii . 

Since N M is related to the residence time by r res = N M / (2dD), we can write 

uj m ~ 2dFL d pDr res = 2dDL d p 2 = uj mf . 



(26) 
(27) 



(28) 



In this way we have shown that for the nucleation rate the mean-field treatment is equivalent to considering particles 
as noninteracting, i. e. counting also meeting events following the first one. For this reason the mean-field value is 
an overestimate of the correct nucleation rate. Furthermore we have proven that 



Wni 

w 



N dSB 



= Af 



(29) 



In Fig. U, the comparison of ui uf /uj, computed exactly in the companion paper pi[ with N aU /N diB , evaluated 
numerically shows clearly that Eq. (^9|) is valid with great accuracy. 

The correction factor Af depends on well known properties of single particles performing a random walk. The 
numerator JV all is just (see Sec. |J): -^.u = 2dL(/3L + a£ ES ). The value of the denominator N diB is well known [^|| in 
the absence of step-edge barriers, being of order L in d = 1 and of order L 2 / InL in d = 2, and it is trivial in the limit 
of infinite barriers, being exactly equal to L d . Hence in d = 1 we obtain Af ~ (L + crf ES ), for all ^ ES . 

In d = 2 we have the limiting expressions TV ~ InL for weak barriers and Af ~ £ E s/L for strong ones. For 
intermediate barriers it is possible to find a simple approximate expression for N dis and therefore an interpolation 
between the two limits. The atom performs on average a number N tr — T res /T tr of traversals of the island. During a 
single traversal each site has a probability pi ~ (1/ InL) to be visited. After all N tr traversals the probability p s that 
a generic site has been visited at least once is given by 1 — p s = (1 — pi) Ntr . Hence we can estimate the number of 
distinct visited sites as 



L 2 



Ps 



L 



l-(l-pi) 



,/Ttr 



(30) 



This expression assumes all traversals to be independent, which is clearly not strictly correct. However it gives the 
right values in the limits £ ES = and ^ ES = oo and for intermediate barriers its accuracy can be tested numerically. In 
Fig. H we have plotted the ratio N all /N diB as a function of ^ ES , for L = 20. The picture shows a reasonable agreement 
between the analytical estimate and the numerical simulation. 



B. The spatial distribution 

In the previous subsection we have shown that MF theory overestimates the nucleation rate by the quantity Af 
because it counts all meetings between two noninteracting adatoms. We are now going to prove that the identification 
of the mean-field approach as a theory of noninteracting particles is valid for the spatial distribution of nucleation 
events as well. We adopt a continuum notation so that a single proof is sufficient to demonstrate that P NI (x) and 
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FIG. 3. Log-log plot of N aU /N dis versus £ ES for d = 2 and L = 20. 



p 2 (x) are proportional in any dimension, for any value of the ES barrier and for any terrace shape. In the regimes of 
strong and infinite ES barrier this result is trivial because both p 2 (x) and P NI (x) are constant. 

We face the problem of determining the quantity P NI (x) for a pair of adatoms, one with initial distribution 
p( 2 \x,0) = p v (x) and the other with the effective distribution p^ (x, 0) = p c!t (x) — ^2 t p^ 2 \x,t) — p(x). 

We can consider the coordinates (xi,x 2 ) of the two atoms (xj. 2 are vectors in a rf-dimensional space) as defining 
the position x — (x\, X2) of a single particle in a space of dimensionality dl — 2d. This particle moves according to the 
diffusion equation dtp = y V 2 p. The factor | appears because of the different time step At employed in describing a 
single walker (At = \/2dD) or two walkers (At — l/2d'D) on a terrace. 

Integrating in time and defining P(x) = dtp(x, t) one finds 



fV 2 P(x) = -p(x,0) . 



(31) 



An interchange of the two particles [p^(x, 0) = p v (x) and p^(x, 0) = p(x)] is equally legitimate and it is useful to 
use a simmetrized form for p(x, 0): 



p(x,0) = 5 [p(x 1 )p u (x 2 ) +p u (x 1 )p(x 2 )} ■ 
Notice that p(x) is also the solution of the equation DW 2 p(x) = —p"(x). Therefore 

f V 2 P(x) - f [p( Xl )W 2 2 p(x 2 ) + p(x 2 )W 2 p(x 1 )} , 
where V| acts on Xi only and V 2 = V 2 + V 2 ,. Hence 

V 2 P(x)=V 2 [p(x 1 )p(x 2 )} , 



(32) 



(33) 



(34) 



i. e. the function x( x ) = -P( x ) — p(x\)p(x 2 ) is harmonic. It is easy to show that x( x ) must be identically zero. 
Hence: 



P(x) = p(x 1 )p(x 2 ) . (35) 

If we set x\ — x 2 = x, the left hand side is just the nucleation probability P NI (x) at point x between two noninter- 
acting adatoms, and the right hand side is the mean-field prediction. Notice that we have not used the explicit form 
of p v (x). Hence, the proof holds for any initial spatial distribution p v (x), so that the equivalence between MF theory 
and the noninteracting particles model is true even if atoms are not deposited uniformly. 

So far we have rigorously shown that the MF result for P(n) is not exact. However, one may wonder whether the 
error introduced by taking into account all meeting events following the first one is expected to be large or negligible. 
We address this issue by evaluating the relative weight of successive encounters for noninteracting particles. 

Let us consider noninteracting particles and define pj as the fraction of times the j-th meeting event actually occurs. 
Clearly po = 1 and pj > Pj+i- Let us also define the normalized distribution for the j-th nucleation event P? [I (n). 
Notice that PJ"(n) = P^ N \n), the distribution for interacting particles. 







The total distribution of nucleation sites is simply 



P NI (n)=5>,-Pf (n) (36) 
i=i 

and the quantities W and W m are given by: 

W = m, (37) 

oo 

W N1 =^P NI (n) = ^Mj • (38) 

n j=l 

If we now introduce the normalized distribution of all fictitious nucleations following the first one: 

pm{n) = g£|ff^>) (39) 

we can write 

oo 

P NI (n) = W Pf» + ^M^ NI H = WPW(n) + P m (n)(W m - W) . (40) 
i=2 

For weak barriers, in d = 1, the weight W of the first term is constant, while the second one diverges as L. As a 
consequence, for large L the distribution P NI (n) is dominated by the contribution of the fictitious successive nucleations. 
In d = 2 the first term goes as 1/lnL while the second is constant. Again, for large L, the contribution of first 
nucleation events becomes negligible. For strong and infinite barriers, W = 1 while W m is infinite, so P NI (n) coincides 
with P™(n). 

In all cases the MF expression for the spatial distribution of nucleation sites [P NI (n)] is dominated for large L 
by the contribution of the fictitious nucleations following the first one. A priori there is no reason for supposing 
that the distribution of the j-th nucleation event is equal to the distribution of the first one, so we expect that the 
difference between the MF spatial distribution and the exact result persists for all values of L. This will be checked 
and confirmed in the following paper pifl. 



V. CONCLUSIONS 



This paper has been devoted to an accurate investigation of the mean-field approach to the problem of irreversible 
nucleation. The main outcome is the proof that MFT is equivalent to a model where particles do not interact and all 
their meetings are counted as fictitious nucleations. 

In the regime of infinite ES barriers, MFT simply breaks down because it predicts a diverging nucleation rate, in 
contrast to the correct value ui — FL d . In the other, physically more interesting, regimes the equivalence of MFT 
with the model of noninteracting particles implies that lo mf overestimates the correct nucleation rate by the factor 
Af = N all /N diB . This ratio has a clear meaning: a diffusing adatom passes Af times on a visited site. It depends on 
single-particle quantities (iV aU , N diB ) whose expressions are well known from the theory of random walks. 

In Table | we summarize the value of the correction factor Af in regimes i and ii and we report the approximate 
expressions for the nucleation rate u>. They are approximate in the sense that numerical prefactors are neglected, but 
they scale correctly with L,D,D',ao. The lattice constant ao has been reintroduced in order to give dimensionally 
correct espressions. Also, we have made explicit the dependence of ^ ES on D and D 1 , so that only basic quantities 
appear. 

The expression lu ~ F 2 L 5 ao/D', valid in two dimensions for strong step-edge barriers has already been given in 
Ref. [fU^ . It is worth repeating that the nucleation rate in this limit does not depend on the diffusion constant D so 
that the nucleation rate can not be promoted by using surfactants. 

Application of MFT is acceptable only in the regime of vanishing barriers in two dimensions, because in this case 
the correction factor [Af — ln(L/ao)] is a small number, for realistic terrace sizes. 

In order to obtain exact expressions for lo it is necessary to have an accurate estimate of or equivalently of N dis . 
W is a function of the terrace size L and of the ES length £ ES : for strong barriers W = 1, while for weak barriers 
W ~ 1 in d = 1 and W ~ 1/ InL in d = 2. So, for realistic values of L, W depends on L and £ ES much more weakly 
than the other quantities appearing inw = FL 2d pW . However, its dependence is not fully negligible: Figs. Ilia and b 
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TABLE I. We report the nucleation rate ui and the correction factor N in the two regimes of weak (i) and strong (u) 
step-edge barriers and in one and two dimensions. The conditions defining regimes i (r re s ^ n r ) and m (rt r -C r re s <C rd ep ) are 
written in terms of the basic quantities L, D,D' , a . 



show that (in d = 2) for L = 20, 14 7 varies by a factor 10 by changing £ ES from zero to infinity and for £ ES = 0, 
W varies by a factor 3 by changing the terrace size from L — 4 to L = 32 (In 32/ In 4 = 2.5). For comparison, the 
quantity p, which appears along with W in the expression for ui, varies by a factor 50 by changing £ ES /L from to 
6 and by a factor 64 by changing L from 4 to 32, for £ ES = 0. The problem of the exact determination of W will be 
tackled again in the following paper. 

A last comment on rate equations and mean-field approximation is in order here. According to the former, the 
nucleation rate is written uj I{E = o\Dp 2 and the latter corresponds to taking <j\ as a constant. In general, o\ is 
defined through the relation $1 = Daip, where $1 is the flux of atoms attaching to an adatom. The resulting relation 
Wre = $1/0 is exact, if $1 is evaluated correctly; for example, we can solve the diffusion equation for a single walker 
on the terrace, where a second walker is taken as an absorbing sink. Since we have shown that the nucleation rate 
can be evaluated assuming an atom as immobile, such treatment is essentially correct. In other words, if the capture 
number is not taken as a constant, the expression u) RB = (T\Dp 2 may give correct results, but this method has nothing 
to do with the usual mean-field approach. 

Finally, with regard to the spatial distribution, we have provided a very general demonstration of the equivalence 
between mean- field and noninteracting particles. We have also shown that the difference between P(n) and P NI (n) 
is not an effect of the finite size of the terrace and it remains true for large L. The full computation of the spatial 
distribution of nucleation events requires the solution of the dynamical problem of two interacting atoms diffusing on 
a terrace. This problem will be solved analytically and/or numerically in the following paper. 
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APPENDIX A: SINGLE PARTICLE ON A TERRACE 

In this Appendix we summarize the behavior of a single particle on a terrace for all values of the Ehrlich-Schwoebel 
length £ ES . 

1. The stationary adatom density 

The discrete evolution equation for a particle in a cubic d-dimensional space is 

Pn(* + 1) = ^ (Al) 



5 



where n + S indicates a neighbour of the site n. If we sum over t and define the quantity p n = Y^tLoPn{t) we obtain 



^ Pn+S - 2dp n 



+ 2d Pn (0) = 0. (A2) 



The terms in square brackets give the discrete laplacian of p n ; therefore the sum p n = J2 t p n (t) is simply the 
solution of the stationary diffusion equation in the presence of the flux 2dp n (0) . 
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In d = 1, for constant p n (0), it is possible to find the exact discrete solution for any value of £ ES , once we remark that 
the general solution of the homogeneous equation is p n = cq + c\n and a part icul ar solution of the nonhomogencous 
equation is p n = —C2n 2 [the factor C2 depending on the constant term in Eq. ([A|)]. 

Boundary conditions are po = ap\ and pL+i = o-Pl, where a = 4s/(l+4s) goes from to 1 as the Ehrlich-Schwoebel 
length 4s varies from to 00. 

The explicit expression of p n is 

Pn = t [ijasL + {L + l)n - n 2 } . (A3) 

Its normalized version is 

pl s y^V = < n i Ewm [£EsL + {L + l)n n2] ■ (A4) 

l^n Pn Ies-^ H ' g 

In a continuum formalism, the equation is Dd 2 p + F = and the solution in <i = 1 is 

= ^ [4s£ + La; - x 2 } . (A5) 

In d = 2 the solution of the continuum equation is as easy as in d = 1 if we specialize to a circular terrace. If L 
now denotes the radius, the solution for g 

p(r) = ^[L 2 + 2L£ ES -r 2 } . (A6) 
We finally evaluate the average density p on the terrace: 



p = j^L(L + 64s) d=l (A7) 
p = -^L(L + 4i ES ) d = 2 [circular terrace] (A8) 



2. The dynamical problem in one dimension 

We now summarize the dynamical behavior of a single particle on a one-dimensional terrace. The two-dimensional 
case is treated in the next Subsection. 

The discrete evolution equation for the particle is 

Pn(« + 1) = \[Pn+l{t) + Pn-l{t)] , (A9) 

with the usual boundary conditions, po(t) = ap\(t) and pL + i(t) = api,(t). The solution is found by separating the 
space and time variables, p n (t) — X(n)F{t): 

F(t + l) _ X(n + 1)+X(n-1) 

~^T~ X - 2Xtn) 0<A<1. (A10) 

The temporal part is F(t) — X t F(0). The spatial part has the general form 

X(n) = A sin(n0) + B cos(n<p) , (All) 

which gives A = cos <fi. The boundary conditions determine the values of A, B and </>. 

In particular, by imposing the boundary condition in n = one obtains B = bA with b — asin0/(l — acos(f>). 
Using this relation and imposing the other boundary condition in n = L + 1 one obtains 



tan(^) = - v " ~ lj °7 o . (A12) 

/I 1 " * 



(a 2 — 1) sin< 
(1 + a 2 ) cos (j> — 2a 

This equation has L solutions which wc label as (f>k with k = 1, . . . , L. Then the general solution is 

L 

Pn(t) = J] cos * (&) x k(n) , (A13) 

k=l 
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with 



X k (n) = A k sm((/) k n) + B k cos((/> fe n) 



(A14) 



Given p n (t) one can compute S(t), the probability that an adatom is still on the terrace at time t after deposition 
(survival probability) : 



n=l 

Another important quantity is the residence time, defined as 

oo 

Tres = ^2t[S(t-l)-S(t)] 



(A15) 



(A16) 



because S(t — 1) — S(t) is the probability that the particle stays on the terrace exactly a time t. It is easy to check 
that 



T res =^2S(t) = ^ 



4=0 



n=l 



t=0 



= p« 



Recalling App. A 1, for the initial distribution p„(0) = p^ = l/L we have, for any value of the ES barrier, 

(£ + l)(L + 2) 



G 



(A17) 



(A18) 



In order to pass to a continuous time we have to multiply it by At — 1/2D. For large L, r res — + 2~)- This 

result agrees with the relation p = Fr res [see Eq. ( A7)]. 



Unfortunately it is not possible to solve explicitly Eq. (A12) for generic values of a: we now discuss the two limit 
cases where an explicit solution is possible. 



a. Zero barriers 



For £ ES — (a — 0), the allowed values of <j>k are 



L + l 



and the general solution is 



Pn(t) = ^ A k COS* 



fc=l 



kit 

L + l 



nkir 
L + l 



with 



.4, 



L 



2 \ ^ . / nkir 

^p„(Q )sm ( 



In particular, two forms of p n (0) are most interesting to us. For a uniform distribution p ra (0) = p v n 
coefficients are 



M = A u k = 



2 (L kir \ fkir\ (l kir 

— — — r sin — sin — esc - 

L(L + l) \2L + l) \2 J \2L + l 



(A19) 
(A20) 

(A21) 
1/L the 

(A22) 



For the distribution p n (0) — Pn — 



L(L + l)(L+2) ' 



(L + 1 — n) [see Eq. (A4)], the explicit solution is: 
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A k = Ai 



6 



sin(f) 



L(L + l)2(L + 2) 8 - 3 



2(L + 1) 



Lkir 



2(i + l) 



(A23) 



As shown in App. A 1, is the normalized version of p„ = X)t=oP™( t ) wnere Pu{t) is the solution of the diffusion 
equation with uniform initial condition p^. Writing explicitly the sum we obtain 



^pu(t) = J2J2 A k cosT f ' sin 



T = 0fc=l 

L 



E 



-1 1 — cos 



fc=l 



(*) 



l + i; \L + 1 

nkir 

sin ' 



L + l 



Hence 



A£ cx 



4f 



4f 



1 — cos 



•(ft) 



2 sin 



2(L + 1) 



as can be easily verified by comparing Eq. ( A22| ) with (A23). 

If we sum p n (t) over n [see Eq. (A20)] we obtain the survival probability: 



(A24) 
(A25) 

(A26) 



S(t) = A k COS* 



k=l 



kn 



sm — 



L k-n 



...... . sm — esc 

L + l \2 L + l 2 V2L + 1 



kn 



1 k-n 



The distribution p n (t) is in general a sum of exponential decays 

L 



s—^ I nkn 

Pn(t) = A k sm ( - ; - | exp 

fc=i 



L + l 



i In cos 



kn 

L + l 



(A27) 



(A28) 



It can be considered as a single exponential when the second slowest decaying exponential is negligible. For L ^> 1 
this means 



exp 



Hence, for t ^> r tJ . 



Pn(*) — Ai si 



For the same reason, for t ^> r tr 



L+l 



27ry * 
L~ J 2 



exp 



L 2 

«1 => t > — - ~ T tr • 
27T Z 



7T \ 2 £ "I / n7r 

— - ~ Ai sm I exp 

LJ 2 V L + 1 / 1 \ 



.. . 2(L + 1) / t 
S(t) ~ Ai — exp I 

\ 7"res 



(A29) 



(A30) 



(A31) 



b. Infinite barriers 



For £ ES = oo (a = 1) the allowed values of 4>k are 

kn 



k = — {k = 0, . . . , L - 1) 



and Afe = Bk tan (§r), so that the general solution is 



L-l 



Pn{t) = Y A k cos* ( ^ ) A fc (n) 



(A32) 



(A33) 
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where 



X k (n) 



kir\ . ( nkn 
tan| 2Zj Sm — 



nkir 
~L~ 



The cocfBcients A)~ depend on the initial condition through the relation: 

L 

.v fe 



1 L 



where (Sko is the Kronecker symbol): 



N k 



1 + tan z 



kit 
2L 



(i + M . 



(A34) 



(A35) 



(A36) 



p n (t) is the sum of a constant (the term for k — 0) and exponentially decaying terms (fc > 0). For p n (0) = l/L, 
the only nonvanishing coefficient is Aq = l/L and this implies for all times 



Pn(t) 



1 



In the general case of nonconstant p n (0), the exponential decays are negligible when exp 
to say t > ~ r tr . 



sL) 2 



(A37) 
<C 1, that is 



c. Strong barriers 



Let us consider now the case of finite but large ^ ES (^ ES L). The solution of Eq. (A12), with a — *■ 1 and large but 
fixed L, yields, for the two smallest 4>k 



2(1 - a) 



L 



fa = 1 + 0(1 -a). 



LL 



(A38) 
(A39) 



The slowest decays in the general solution are therefore exp(— </>ii/2) and exp(— <f)\t/2). For finite values of L we 
can neglect the second exponential for times such that 



exp [tin cos (02)] <C 1 
Hence for times larger than r tr one can write 



p n (t) = Bi cos n 



cos 



LL 



B\ exp — 



LL 



(A40) 



(A41) 



where B\ ~ l/L and r res = L£ ES . This value of r res , multiplied by At = 1/2D coincides with its continuum 
counterpart p/F — L^ ES /(2D). 

3. The dynamical problem in two dimensions 

It is useful to summarize here some results for a single particle on a two-dimensional terrace. The general solution 



P„ 



1 



k,j=l 



kn 
L + l 



L+l 



X k {m)X 3 {n) 



(A42) 



where the coefficients Akj are 
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1 

A kj = J^-jy- ^ Vm,n{ Q ) X h{m)Xj{7 



(A43) 



m,n— 1 



For zero barriers -Xfc(n) = sin f f^-J and Nu = For a uniformly distributed adatom, p„ „ = 1/L 2 and the 

coefBcients are 



A u _ A u A u 

A kj — A k A j 



-i 2 



Ljw 



sm 



/i'TT 

T 



2(L + 1) 



CSC 



sin 

A'TT 



J7T 

2 



sm 



2(L + 1) 



CSC 



Lkn 
.2(1+1) 



2(L + 1) 



(A44) 
(A45) 



We indicate as n the normalized solution of the stationary diffusion equation in the presence of a constant flux. 
Differently from what occurs in the one-dimensional case, the explicit form of p^ n is not known exactly for a square 
terrace. However, the expression of its coefficients Afj can be obtained by exploiting the property (see App. A 1 ) 
that pfn n — N^2^°_ p mj7l (T) where Pm,n(i~) is the solution of the diffusion equation with uniform initial condition 
Pm n = and ./V = l/r res is a normalization factor: 



oo oo L ^ 



T = 



t=0 k,j=l 
L 



k~K 

L + l 



L+l 



feli= il-i [cos(^)+cos(^j 



mkir 
L + l 



mkir 
L + l 

njir 



L + l 



L + l 



(A46) 



Hence 



A * = N - 



A kj 



cos( 7 ^ r )+cos(^ T ) 



(A47) 



The numerical prefactor N can be determined by imposing that the sum over m and n of p^ n is 1, that is 



X] X! A kj 

m,n k,j 



sm 



mkir 
L + l 



sm 



njir 
L + l 



= 1 



(A48) 



which implies 



1 

N 



L{L + 1) 



k,j L 2 



( A l? 



cos 



f J=E_^ + COS (j^- 



U+i 



U+l 



(A49) 



In the limit of large L, 



N 



L(L + 1) 



1 2 



1 — COS 



2 7 2 



(A50) 



Hence in the continuum r, 



, 32 l - 
res — AlKl 



In the limit of strong but finite barriers one finds r res = L£ ES /(4D), so that a = 1/4. 



(A51) 
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APPENDIX B: THE EFFECTIVE DISTRIBUTION 



We want to evaluate the effective distribution: 

oo 

^E P -p(tK(t) 
r=0 



(Bl) 



introduced in Sec. Ill A. Since Pu{t) decays to zero after a time of order r res , for regimes i and ii (where r res <C Tdep) 
Pde P (i~) can be taken as a constant. Hence 



o°n = *^2Pn( T ) ■ (i and ii 

T dep T=Q 



(B2) 



The sum ^ r p n (r) has been shown in App. A 1 to be equal to the solution p n of the stationary diffusion equation, 
which always has a parabolic shape. Its normalized version is 



pn 



Pn 



En=l Pn Tres 



so that 



r ff 'res s 

Pn = —Pn 

Tdep 



(i and ii) 



(B3) 



(B4) 



This equation corresponds to Eq. ( |l4| ) in the limit r res <C Tdep- 

For strong and infinite barriers (regimes ii and Hi), the contribution of times shorter than r tr is smaller than T tr /Tdep 
and therefore negligibl e. Ac cordingly, we can evaluate the sum (Bl) using the expression for p„(r) that is valid in the 
limit r ^> T tr [see Eq. (A.41), in d — 2 the generalization is triviaTf7p n (''") = (l/^ d ) exp(— r/r res ) and obtain 



Pn 



I'' 



7~ H exp 



T dep T=Q 

Converting the sum over discrete times into an integral we have 



(- + - 

\ Tdep <res 



1 



Pn 



-Pn > 



Tres ~h Tdep L T res -f- Tdep 

where - as usual - p^ is the normalized solution of the stationary diffusion equation. Thus, formula (Jl4|): 



Pn 



' res s 



Tres T" Tdep 



(B5) 



(B6) 



(B7) 



can be used in all the different regimes. 

In the continuum it is possible to work out a more rigorous approach and determine pjf as the solution of a 
differential equation, which is the generalization of the stationary diffusion equation (^). We start with the diffusion 
equation for p(x, t), dtp = D\7 2 p, where a; is a d— dimensional vector. If we multiply both sides by P dcp (t) and integrate 
in time, we obtain: 



/>oo />oo 

/ dtP dep (t)d t p(x,t) = DV 2 dtP dep (t)p(x,t) 
Jo Jo 

The right hand side is just DV 2 p Bll (x) while the left hand side is: 

dtP dep (t)d tP (x,t) 

dt[d t P dep (t)]p{x,t) 
— f dtP dep {t)p(x,t) 

Tdep Tdep JQ 



P dQp (t)p(x,t)\™ 

1 1 



(B8) 

(B9) 
(BIO) 
(Bll) 
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and Eq. (fBq) becomes: 



HVVfi) P c!t (x) + F = 

Tdep 



(B12) 



It differs from the stationary diffusion equation DV 2 p + F = because of the presence of a 'desorption' term 
( — P clt (x) / Tdep) which is the responsible of the saturation of p° li (x) at large £ ES . As a matter of fact, in the limit 
£ ES — > oo, p is known to diverge as j^L£ ES [see Eqs. (A7,A8)] while the above equation clearly shows that p' !t (x) goes 
to the constant FTd ep = 1 A4. 

The exact solution of ( |B12 ) can be found both in d = 1 and in d = 2 for a circular terrace and the proof that 
p Bti (x) = — — p s (x) works much in the same way in the two cases. We give here a few more details for the 

bidimensional case. The solution of Eq. ( B12| ) with the usual boundary condition d r p Bti (r) = —p° !! (r) / '£ ES evaluated 
for r = L (the radius of the circular terrace) is 



P° i! (r) = -4a 



1 - 



y/rTFTl 



/o(- 



-) + 



-h(- 



■) 



(B13) 



' \J DT dcp y/ Dr dep y/ Dr dc 

where Iq and I± are the modified Bcssel functions of order zero and one, respectively. The arguments of the Bessel 



functions are at most equal to ( 



y/D^Z 



8r tr 



p^(r) 



a small quantity. An expansion of the Bessel functions gives: 
L 2 + 2£ ES L - r 2 



ttL 2 AD Tdep + L{L + 2£ es ) 



By using the results (|A6|) and (A8), after some algebra we obtain the final expression: 



p rff (r) = 



T dep +£1(4 + 



2 > 



P s (r), d = 2 



The calculation in d = 1 leads to the result 

p^(x) 



Tdep +£>(«+ 2 S ) 



p s (x) 



1 . 



(B14) 



(B15) 



(B16) 



The quantity ■ •) appearing on the right hand side at denominator does not coincide with r res because the term 
has a prefactor j instead that ^ in d — 2 and a prefactor | instead that ^ in d = 1. Nonetheless such quantity 



differs from r r , 



for a quantity of order rj r which can be safely neglected with respect to Tdep (always appearing at 



denominator) so that, in the limit r tr <C Tdep ( a limit applied throughout the paper) we can conclude that the relation 

P act (x) = Tr : s _ P s (x) (B17) 

is always valid. 



Tdep 
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